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Abstract 

We analyze the constraints on neutrino mass spectra with extra sterile neutrinos as 
implied by the LSND experiment. The various mass related observables in neutri¬ 
noless double beta decay, tritium beta decay and cosmology are discussed. Both 
neutrino oscillation results as well as recent cosmological neutrino mass bounds are 
taken into account. We find that some of the allowed mass patterns are severely 
restricted by the current constraints, in particular by the cosmological constraints on 
the total sum of neutrino masses and by the non-maximality of the solar neutrino 
mixing angle. Furthermore, we estimate the form of the four neutrino mass matrices 
and also comment on the situation in scenarios with two additional sterile neutrinos. 
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1 Introduction 


Scenarios with four neutrinos became popular on the wake of the LSND evidence of — v e 
transitions JET]. Interpreted in terms of neutrino oscillations, the indicated mass scale 
for LSND is in the eV 2 range. Together with the evidence for neutrino oscillations from 
atmospheric (plus K2K) and solar (plus KamLAND) neutrino observations, requiring mass 
scales around 10 -3 eV 2 and 10 -4 eV 2 , respectively, a fourth sterile neutrino has to be 
introduced in order to accommodate the presence of three distinct mass squared differences. 
A priori, four neutrino scenarios allow for two possible mass patterns: 

(i) 2+2 scenarios, in which two pairs of neutrino states are separated from each other 
by the LSND mass scale. There are two possibilities for 2+2 scenarios; 

(ii) 3+1 scenarios, in which one single neutrino state is separated by the LSND mass 
scale from the other three states. There are four possibilities for 3+1 scenarios; 

Oscillation analyzes in both schemes were performed by a number of authors ElEUSE] and 
also the astrophysical and cosmological implications were investigated Eld- Historically, 
among the above two alternatives the 3+1 scenarios were at first relatively disfavored ;3J 
because of the non-observation of oscillations in short baseline experiments like KARMEN 
0, Bugey (HI and CDHS pTO]. Therefore the 2+2 scenarios were found to be more com¬ 
patible with the existing data. The sterile neutrino oscillation solution in 2+2 scenarios 
was viable for both solar and atmospheric neutrinos. However, SuperKamiokande data 
disfavored oscillation of the atmospheric v [L to purely sterile neutrinos 0, and later on 
the SNO data started establishing the neutral current component in the solar v e flux P- 
For some time a mixed scenario, where the atmospheric neutrino anomaly is due to z/ M — v S)T 
and the solar neutrino anomaly is due to v e — u SjT , remained compatible with all data fTr>1 . 
However, all recent analyzes show that 2+2 scenarios are ruled out at a high a from the 
existing data mm- Both atmospheric and solar neutrino data strongly disfavor oscil¬ 
lations to pure sterile species. This disfavored the 2+2 scenarios irrespective of whether 
LSND results are confirmed or not. The most updated analysis in the 3+1 scheme per¬ 
formed in mm shows that non-evidence of neutrino oscillation in other short baseline 
(SBL) experiments combined with atmospheric neutrino data from SuperK and K2K is 
inconsistent with the LSND signal at 95% C.L. and only marginal overlaps are found at 
99% C.L. Thus, with increased precision of solar and atmospheric neutrino flux measure¬ 
ments the four neutrino explanation of the LSND anomaly suffered a setback. This led to 
many alternative explanations of the LSND anomaly including introduction of two sterile 
neutrinos - the so-called 3+2 scenario |TB] -, CPT violation [IT], quantum decoherence 
effects violating CPT jTS|, mass varying neutrinos ra, neutrino decay in four neutrino 
scenarios m, lepton number violating muon decay 123 , decay of a heavy neutrino ra or 
extra dimensional aspects 0. 

Oscillation experiments can only measure the mass squared differences but not the absolute 
masses. The most direct and model independent way to measure the absolute masses is 
via kinematic measurements involving nuclear beta decay. The best bound at present is 
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mg < 2.3 eV (95% C.L.) coming from the Mainz tritium beta decay experiment [ 21 ]. The 
KATRIN experiment is expected to increase the sensitivity down to ~ 0.2 eV [22]• 
Information on absolute masses can also come from neutrinoless double beta decay (Ou/3/3). 
Neutrinoless double beta decay experiments aim at observing the process 

(A, Z ) —■> (A, Z + 2) + 2 e 

This is a lepton number violating process and its observation will establish the Majorana 
nature of neutrinos j2E]- The decay width depends quadratically on the so-called effec¬ 
tive mass. We assume here that only the light Majorana neutrinos implied by neutrino 
oscillation experiments are exchanged in the diagram of 0 u/3/3. In the basis in which the 
charged lepton mass matrix is real and diagonal, the effective mass is then nothing but the 
absolute value of the ee element of the neutrino mass matrix. The best current limit on 
the effective mass is given by measurements of ' 6 Ge established by the Heidelberg-Moscow 
collaboration [2Zj (with similar results obtained by the IGEX experiment 2BJ) 

(m) < 0.35 ( eV , (1) 

where ( = 0(1) indicates that there is an uncertainty stemming from the nuclear physics 
involved in calculating the decay width of 0 u/3/3. The running projects NEM03 j22] and 
CUORICINO [SU] will be joined in the near future by next generation experiments such 
as CUORE EH, MAJORANA EH, GERDA EH, EXO EH, MOON EH, COBRA EH, 
XMASS, DCBA [27;, CANDLES EH, CAMEO -{0 (for a review see [40])• One can safely 
expect that values of (m) one order of magnitude below the limit from Eq. © will be 
probed within the next, say, 10 years 1 . This means that scales of order a/A mLg ND will be 
fully probed, and are even under investigation now. Since the effective mass measured in 
0 u/3/3 also depends on the neutrino mixing angles, the neutrino mass scale and ordering, 
as well as the mass squared differences, it is possible to obtain additional constraints on 
sterile neutrino scenarios using neutrinoless double beta decay E21 EH EH- 
Important constraints on sterile neutrinos can also come from cosmology. Inclusion of 
an extra neutrino, even if sterile, can be in conflict with cosmological observations. The 
problems are increased if the extra sterile neutrino is massive and has significant mixing 
with the active species. In particular the Big Bang Nucleosynthesis model of standard 
cosmology, which explains light element abundances of the Lhiiverse, puts constraints on the 
number of neutrino species. The latest bound found in [45; for instance is 1.7 < N v < 3.0 
at 95% C.L. and in ESj A is quoted that N v = 3.147o;g5- The differences in the results 
are due to different inputs regarding the uncertainties in the primordial He abundance. 
Observations of the Cosmic Microwave Background and of large scale structures can also 
constrain the number of neutrino species. A summary of these bounds obtained by various 
groups including different data sets can be found in A7J. The upper limit on the number 
of neutrinos in these analysis can vary from 6 to 8. A recent bound as quoted in [47J is 

1 Not to forget, those experiments aim also to put the controversial m evidence of part of the 
Heidelberg-Moscow collaboration to the test. 
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N u = 4.2+1 2 a f 95% C.L. Another important constraint from cosmology comes on the sum 
of total masses of all the neutrinos, £ = For four light neutrinos with degenerate 

masses the bound is £ < 1.7 eV (95% C.L.) from WMAP and 2dF data |i7j . For four (five) 
neutrinos, with one (two) of them carrying a mass, the bound is £ < 1.05 (1.64) eV (95% 
C.L.) HE]. Improvement of these numbers within one order of magnitude is expected 03- 
Note that these bounds depend on the priors and data sets used, for slightly more stringent 
bounds see, e.g., [32] • The above constraints can however be evaded if the abundances of 
sterile neutrinos in the early Universe can be suppressed. This requires going beyond the 
framework of standard cosmology and introducing mechanisms such as primordial lepton 
asymmetries ra, low re-heating temperature 03. additional neutrino interactions jE2J 2 
etc. 

Turning back to oscillations, the MiniBooNE experiment |3) is expected to confirm or 
refute the LSND signal and is expected to publish results within the next 6 months or so. 
If MiniBooNE does not confirm the LSND signal, then with the data collected with 10 21 
protons on target they can rule out the entire 90% area allowed by LSND with 4 to 5a 
. If however they confirm the LSND signal then this will give rise to an intriguing sit¬ 
uation in what regards the explanation of global oscillation data from accelerator, reactor, 
atmospheric and solar neutrino experiments. If confirming the LSND result, MiniBooNE 
can not distinguish the allowed four (or Eve) neutrino mass spectra. To understand the 
implied mass and mixing scheme, other observables are therefore crucial. This concerns in 
particular observables depending on the neutrino masses and ordering. Inasmuch one can 
use these future measurements to identify the neutrino spectrum is one of the motivations 
of this work. We stress here that we assume only the neutrino oscillation explanation 
of the LSND result is correct, i.e., the new physics alternatives (not necessarily predict¬ 
ing a signal for MiniBooNE) put forward in Refs. [T71IT%irH?ll2Ull2Tll221l23| are not required. 

In this paper we examine what constraints from current and future data can be obtained 
on possible neutrino mass spectra in scenarios with one or more sterile neutrinos. For the 
four allowed 3+1 scenarios we give the neutrino masses, their sum as testable in cosmology, 
the kinematic neutrino mass for tritium experiments and the effective mass in neutrinoless 
double beta decay. We include the most recent values of mass-squared differences and 
mixing angles from latest global analyzes of oscillation data. We furthermore reconstruct 
the possible mass matrices in four neutrino scenarios that are consistent with the current 
data. Finally, we also comment on 3+2 scenarios. 

The paper is build up as follows: In Section |2] we discuss our parametrization of the four 
neutrino mixing matrix and summarize the relevant formulae for the neutrino masses, their 
sum, the kinematic neutrino mass measured in beta decay experiments and the effective 
mass that can be observed in neutrinoless double beta decay. In Section 01 we apply this 
framework to 3+1 scenarios. Approximate forms of the neutrino mass matrices in 3+1 
schemes that are consistent with the current data are given in Section 0 In Section 01 
we comment on the above quantities in the 3+2 scheme, before presenting our summary 

2 See however G31 
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and conclusions in Section El The oscillation probabilities for the relevant short baseline 
oscillation experiments are delegated to the Appendix. Although 2+2 scenarios are highly 
disfavored we also add for the sake of completeness an Appendix on the implications of 
such scenarios for neutrino masses from cosmology, beta decay and neutrinoless double 
beta decay. We also discuss the form of mass matrices in the 2+2 scenarios. 


2 Four Neutrino Mixing and Neutrino Masses 


Neutrino mixing is described by the unitary Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 
matrix U $ 2 ]. Fc> r four Dirac neutrinos it contains 6 angles #12,13,14,23,24,34 and three 
phases £ 13 , 14,24 (“Dirac phases”). Their Majorana nature, which we shall assume, adds 
another three phases (“Majorana phases”), which do not have any consequences in neutrino 
oscillations (SH)- We parametrize U as 
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(3) 


with the usual notation s tJ = sin 0 i3 and = cos 0 i3 . The diagonal matrix P contains the 
three Majorana phases, which we denote a, (3 and 7 : 


P = diag (1, 


a -ia/2 -i(f3/2-6 13 ) 


3—*(7/2—+4 A 


(4) 


For most purposes it is sufficient to analyze the individual experimental data in a two-flavor 
framework. Depending on the neutrino mass spectrum, one can then identify certain ele¬ 
ments of the PMNS matrix with the mixing angle in a two-neutrino oscillation probability. 
For the parameters governing solar (and KamLAND), atmospheric (and K2K) and short 
baseline reactor neutrino oscillation it holds at 3cr EH EE] 

7.0 x 10~ 5 eV 2 < A m% < 9.3 x 10" 5 eV 2 , (5) 

0.25 < sin 2 0 Q < 0.40 , (6) 

with best-fit values of A m 2 Q = 8 x 10 -5 eV 2 and sin 2 0 & = 0.31. The atmospheric mass 
squared difference and mixing angle at 3cr are known within ra 

1.3 x 10 -3 eV 2 < A m 2 A < 4.2 x 10 -3 eV 2 
0.33 < sin 2 6 *a < 0.66, 
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Figure 1: The four allowed 3+1 mass orderings. 
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with best-fit values of A m\ = 2.2 x 10 3 eV 2 and sin 2 6 *a = 1/2. The mixing angle $13 at 
3a is restricted to lie below the value H 

sin 2 ^chooz < 0.044 . (9) 

In Appendix 0 we give the expressions for the oscillation probabilities in short baseline 
accelerator and reactor experiments and also for the 1 km reactor experiment CHOOZ. 
For all the short baseline experiments excepting CHOOZ the one mass scale dominance 
approximation holds to a good precision. A comparison of the probabilities in the 3+1 
picture (recall that 2+2 scenarios are highly disfavored) with the two generation LSND 
probability reveals that sin 2 26* LSND in 3+1 scenarios is always of the form 4\U ei \ 2 \ U^ 2 , 
where U ei | 2 is constrained to be small from Bugey reactor and solar neutrino data and 
| | 2 is constrained to be small from CDHS and atmospheric data. The index i depends 
on the mass ordering. In Ref. nmg the allowed area in the 3+1 schemes is given in the 
Am 2 SN D — sin 2 2#lsnd plane. The plot shows two overlap points at 99% C.L. corresponding 
to 


(Am£ SND , sin 2 20 LSN d) = (0.9 eV 2 , 0.002) and (1.8 eV 2 , 0.001) . (10) 


6 


























LSND 


LSND 



2+2A 2+2B 

Figure 2: The two allowed 2+2 mass orderings. 

In our following analysis we take these two allowed values as the illustrative values of 
AmLSND and take \U e i \ 2 as 0.008 and 0.004, respectively 3 . It is to be noted that the 
MiniBooNE sensitivity plots in, e.g., Ref. give an allowed region around AmL SND = 0.9 
eV 2 , which would be obtained if they confirm the LSND signal 0 - Another allowed region 
is found for which 

A m 2 LSND = (0.2 - 0.5) eV 2 and sin 2 20 LSN d = 0.01 - 0.04 . (11) 

We therefore give in the following another set of plots for which we allow Am 2 SND and 
sin 2 26*lsnd to vary in this range. To extract \U e i \ 2 from sin 2 2$lsnd we assume that 
| Dei | ~ U /n | . With this approximation we have 4 |D ei | 4 = sin 2 26 ) LSND, which gives \U e i \ 2 ~ 
0.05 — 0.1. It is to be noted that this range is not allowed at 99% C.L. according to the 
analysis of m- Nevertheless we take this range as an illustrative example to compare with 
the other two cases with a relatively higher Am 2 SND and lower sin 2 26 l L g ND . 

The neutrino mass matrix is given by 

m v = U* m* ag I/t ) where m)) ias = diag(mi, m 2 , m 3 , m 4 ) . (12) 

We will order the four neutrinos such that m .4 > m 3 > m 2 > mi so that Am 2 SND is always 
given by Am 2 SND = m\ — m\ = Am^. Depending on the relative ordering within the 

3 This corresponds to |R P i| 2 = 0.065, which is the highest allowed value of at 99% C.L. according 

to un. 
















Figure 3: The individual neutrino masses 777 .^ 2 , 3,4 and their sum £ as a function of the 
smallest neutrino mass m We have 3+lAa and 3+lAb on the left (they generate basically 
identical results), 3+1B in the middle and 3+1C on the right. We chose Am 2 SND = 0.9 
eV 2 and fixed A m% and A m\ to their best-fit values. 


scheme, there are four possibilities for 3+1 and two possibilities for 2+2. We display all 
six cases in Figs. |^and[2l One might compare this situation with the three-flavor case, 
where there are only two possibilities, the normal and inverted ordering. The other two 
mass differences correspond to A/77 q and A m\. One can determine the individual masses 
as functions of the smallest mass mi and the three mass squared differences: 

m 2 = y/m 2 + A 7772 ! > 

7773 = a/ m\ + Am^ + Am| 2 , (13) 

7774 = y/mf~+ A rnli + Am§ 2 + Am | 3 = y/mf + Am 2 SND . 

As an immediate application, we can then calculate the sum of neutrino masses £, 

S = mi + 777 2 + m 3 + m 4 . (14) 

for which interesting constraints from cosmology apply. Since the individual neutrino 
masses mi, 2 , 3 , 4 depend crucially on the mass spectrum, their sum £ will do so as well. For 
the four possible 3+1 neutrino spectra from Fig. |T| we display the neutrino masses and 
their sum £ as a function of the smallest mass in Fig. [j] 
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When in addition the mixing matrix elements of the PMNS matrix are specified, one can 
determine mg , the parameter measured in the direct neutrino mass searches in nuclear 
beta decay experiments such as KATRIN. We will denote this parameter the “kinematic 
mass”. It is given by: 

mg = \j\U el \ 2 m\ + \U e2 \ 2 m 2 2 + \U e3 \ 2 m\ + \U e4 \ 2 m\ . (15) 

With the inclusion of mixing, the parameters entering the mass measured in beta decay 
can be expressed in terms of the lowest mass, the mixing matrix elements \U e i \ 2 and the 
mass squared differences. Consequently this quantity can also put constraints on the 
possible mass schemes and their ordering [5Dj. For three neutrino frameworks this has 
no observable effect since the future sensitivity on mg corresponds to quasi-degenerate 
neutrinos, for which unitarity of the PMNS matrix leads to no dependence on the mixing 
matrix elements and for which the normal and inverted ordering generate identical results. 
In the various four neutrino scenarios to be discussed in the following, this will change. 
Neutrinoless double beta decay experiments are sensitive to the effective mass which is 
given as 

(m) = = \\U el \ 2 mi + \U e2 \ 2 m 2 e ia + \U e3 \ 2 m 3 e l/3 + |£/ e4 | 2 m 3 e* 7 | 

= | mi Ci2 C13 cf 4 + m 2 e ia cj 3 c \ 4 + + m 3 e z/3 c \ 4 s 43 + m 4 e* 7 s \ 4 | ( 16 ) 

= cj 4 \cj 3 (m 1 cj 2 + m 2 e ia sj 2 + m 3 e l(S tj 3 ) + m 4 e* 7 tj 4 1 . 

We defined here tij = tan Q ir As can be seen, (m) is sensitive to the Majorana phases 
which may be present in the neutrino mass matrix. The three Dirac phases do not appear 
in (m). The effective mass depends on 10 out of the 16 parameters of the general 4x4 
neutrino mass matrix. This might be compared with the three-flavor case, in which (m) 
depends on 7 out of a total of 9 parameters. Moreover, as in the three-flavor case (for 
recent analyzes, see EDE2I), there is a strong dependence on the mass spectrum. 

We conclude that the three mass related observables (m), mg and £ are powerful tools to 
discriminate among the various possible mass orderings. This will be the subject of the 
next Section. 


3 Neutrino Masses and Neutrinoless Double Beta De¬ 
cay in 3+1 Scenarios 

In the next Subsections we discuss the predictions for the sum of neutrino masses, the 
neutrino mass measured in nuclear beta decay experiments and the effective mass measured 
in neutrinoless double beta decay in the different 3+1 scenarios. A common feature of the 
effective mass is that it can be expressed as a known three-flavor contribution obtained, 
e.g., in m, mi plus an additional term related to the LSND scale. 
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3.1 Neutrino Masses and Neutrinoless Double Beta Decay in 
Scenarios 3+lAa and 3+lAb 


As can be seen in Fig. [3 the scenarios 3+lAa and 3+lAb have three quasi-degenerate 
neutrinos with a mass given by the LSND scale and a fourth, lightest state separated by 
the LSND scale. It holds that Am|, = mf — m A m\ = m| — m 2 and Am = m\ — m 
A m\ = m'l — m§, respectively. For the 3+lAa case we can use Eq. m to express the 
masses in terms of the smallest mass and the three mass squared differences as 

m 2 = sjm\ + Am£ SND - Am| - A m\ , 

m 3 = y/m{ + Am£ SND - A m 2 & , (17) 

m 4 = y/m\ + Am^ SND . 


For the 3+lAb case Am^ and Am 2 ., replace each other. Since the neutrino masses are 
governed mainly by Am 2 SND , the predictions of scenarios 3+lAa and 3+lAb for all mass 
related observables are almost identical M and therefore we treat these two cases together. 
Since the minimal mass of the three quasi-degenerate neutrinos is a/A mL SND , it follows 
for the sum of neutrino masses that 

£3+lAa,b > 3 -^/A'/TT. 2 gND ■ (18) 

The left panel of Figure 01 shows the four masses and their sum £ as a function of the 
smallest mass mi. The three heavier masses are indistinguishable and much larger than 
the lightest state unless mi ~ 1 eV. In Figure 01 we have fixed Am 2 SND at 0.9 eV 2 and 
hence y 3 + 1Aa > b > 2.8 eV. It follows that these two schemes would already be in conflict 
with a cosmological limit of £ = 1 eV unless Am 2 SND ^0.1 eV 2 . Nevertheless, let us 
discuss these scenarios further. One would expect that mi -C m 2 , 3,4 holds, so that the 
fourth state with mass mi effectively decouples in what regards the predictions. 

Neutrino data implies in scenario 3+lAa (3+lAb) that sin 2 20 lsnd = 4 \U e i\ 2 |C/ M i| 2 (see 
Appendix 1X1) and that sin 2 ^chooz — |L/ e2 1 2 (sin 2 ^chooz — |C/ e4 | 2 ). One furthermore has 
\U e3 \ 2 ~ cos 2 9q and |t/ e4 | 2 ~ sin 2 0 e (case 3+lAa) or \U e3 \ 2 ~ cos 2 d Q and \U e - 2 \ 2 — sin 2 0 Q 
(case 3+lAb). 

To an excellent approximation, the predictions for neutrinoless double beta decay are very 
similar to the three-flavor case with quasi-degenerate neutrinos |44j . In that case the normal 
and inverted mass ordering can not be distinguished via Ov/3/3, since their predictions differ 
only by corrections of order Am^/mo, where mo is the common three-neutrino mass scale. 
Analogously, the scenarios 3+lAa and 3+lAb can not be distinguished via Ou/3/3, since 
they generate identical results up to corrections of order A'm)(/Am 2 SND . For 3+lAa, one 
has 
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(19) 
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Figure 4: The effective mass as a function of the lowest mass in scenario 3+1A. The left 
most column is for A?n 2 SND = 0.9 eV 2 and the middle column is for Am 2 SND = 1.8 eV 2 . 
The right most column is for Am 2 SND varied between (0.2 - 0.5) eV 2 . The other parameters 
are varied in their current 3a allowed range and all the phases are varied between 0 and 27 t. 
Also shown is the mass mp that will be measured in beta decay experiments, the current 
and a prospective future limit on the effective mass and the future KATRIN limit. 


We have defined c Q = cos# Q , s Q = sin 9 Q and used that Am 2 SND A m\ Am 0 . We 
also have defined 


(m)^ D = m 0 (cos 2 9 q + e^ Q ^ sin 2 9 Q + e ^ 7 ^ sin 2 ^chooz) 

with |(m)? D | < mo 1 - ta n 2 g 12 ~ 2 siir flcuooz ^ mQ cos 2 fl Q ? 

1 + tan 9 q 


( 20 ) 


with m 0 = sj AmL SND + m\ . The quantity ( 771 ) 3 : ° J us ^ the usua l three-flavor effective 
mass with the common mass scale given by the LSND scale. The contribution of (m)^ D 
is at least one order of magnitude above the term corresponding to the LSND scale: for 
instance, if mi = 0.3 eV, Am 2 SND = 0.9 eV 2 , and with cos2 # 0 ~ 0.2 — 0.5, one has 
(m)? D 0.2 eV and rri\ sin 2 2d LSND ~ 0.002 eV. 
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Hence, neglecting further m\ with respect to Am 2 SND , and also sin 2 6 l cHooz, one finds 

(m) 3+1Aa ~ yj A ml SND \J 1 - sin 2 29 Q sin 2 (0 - 7)/2 . ( 21 ) 

Case 3+lAb is obtained by replacing a with 7 . Due to the non-maximal solar neutrino 
mixing, the effective mass can not vanish in case of scenarios 3+lAa and 3+lAb, its range 
is given by \J Am 2 SND cos20 o ^ (m) 3+1Aa,b ^ \J Am 2 SND . Therefore, for Am 2 SND = 0.9 
eV 2 it follows that (m) min ~ 0.19 eV whereas for Am 2 SND = 1.8 eV 2 we have {m) min ~ 0.27 
eV. If AmL SND is in the range (0.2 - 0.5) eV 2 , then (m) min will be somewhat lower. This 
approximate behavior is reproduced in Figure 0] where we have plotted the effective mass 
as a function of the smallest mass. Also included in the Figure is the predicted value 
of the parameter mp , together with the anticipated KATRIN limit on mp of 0.2 eV. For 
comparison, we also gave the present Heidelberg-Moscow bound on the effective mass, 
together with a prospective future limit of 0.04 eV. Note that the minimum values of (m) 
in the limit of vanishingly small m\ are slightly lower than those mentioned above because 
of the non-zero value of sin 2 #chooz- The limit on the effective mass of 0.35 ( eV rules out 
part of the predicted range: for instance, if sin 2 d Q = 0.28 and Am( SND = 0.9 eV 2 , then 
sin 2 (0 —y)/2 > 0.93 to obey the constraint of (m) ^ 0.35 eV. For sin 2 # Q = 0.31 (0.25; 0.40) 
and Am 2 SND = 0.9 eV 2 , we have that (m) ^ 0.36 (0.48; 0.19) eV, whereas for Am 2 SND = 1.8 
eV 2 it holds (m) ^ 0.51 (0.69; 0.27) eV. Therefore, if sin 2 0 Q turns out to be on the lower 
side of its currently allowed range, then scenarios 3+lAa,b face serious problems with the 
constraints from 0 u/3/3. 

The kinematic neutrino mass in scenarios 3+lAa,b is directly given by the LSND scale: 

m 3+lAa,b _ A m 2 gND _ (22) 

For fixed Am 2 SND the prediction for mp is therefore a line, whereas when a range of values 
for AmL SND is given, we also have a range of values for mp. 

We can summarize the situation 1 for scenarios 3+lAa,b as follows: 

(m) 3+1Aa,b ~ m^ +1Aa,b y 1 — sin 2 2 9 Q sin 2 0 ~--- \J 1 — sin 2 2 d Q sin 2 0 , (23) 

where 0 is some combination of Majorana phases. Some constraints on these parameters 
might be obtained in this scenario. The kinematic mass mp predicted by this scenario 
is much above the KATRIN sensitivity and it will be disfavored if KATRIN confirms mp 
around 0.2 eV. The same is true when cosmological searches do not find a signal close to 
their current bounds. Future limits on the effective mass below roughly 0.05 eV will also 
rule out scenarios 3+lAa,b. A limit of 0.1 eV will rule out the two overlap points at 0.9 
and 1.8 eV 2 . 

4 We note here that the predictions for Ou/3/3 and mp are identical (i.e., up to corrections of order 
A?71 a/Atolsnd) to the one of (the highly disfavored) scenario 2+2A, treated in Appendix El Therefore, 
one can not distinguish these scenarios via (D/3/3 or tritium decay experiments. However, since one has 
£ 2 + 2 A ^ 2 £3+iAa,Y there would be a chance to distinguish them via cosmological measurements. 
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Figure 5: Same as previous Figure for scenario 3+1B. 

3.2 Neutrino Masses and Neutrinoless Double Beta Decay in 
Scenario 3+1B 

The structure of scenario 3+1B is depicted in Fig. |T| We identify A uIq = — m\ and 

A m\ = m\ — m\. The heaviest neutrino 7774 is separated by the LSND scale from the 
remaining three, which enjoy a normal ordering. We can express m 2 , m 3 and 7774 in terms 
of the lowest mass rri\ and the three mass squared differences as 

777-2 = \/m\ + A'/77q , 

7773 = \Jm\-Y Am 2 A + A , ( 24 ) 

7774 = sj 777^ + A?77^ SND . 

In Fig. 01 we show - with A?77 2 snd taken as 0.9 eV 2 - the four masses as well as their sum as 
a function of the smallest mass 7774 . We have “unification” of tt 7 2 and 7774 when 7774 ^ 0.01 
eV and of 7773 and ??7 2 when 7774 ^ 0.1 eV. For small 7774 ^ 0.01 eV one finds 

y 3 + 1B ~ ^A/77 2 snd . (25) 

To have S < 1 eV, it is required that 7774 < 0.01 (0.1) eV if A/ 77 2 SND ~ 1 (0.2) eV 2 . 
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In this scheme it turns out that \U e ^\ 2 ~ sin 2 0 chooz and 4 |[/ e4 | 2 \U^\ 2 = sin 2 20LSND- 
The remaining elements of the e-row of U are \U e i\ 2 ~ cos 2 0 0 and |C/ e 2 1 2 — sin 2 0 0 . Ne¬ 
glecting mi, we have m 2 — m 3 ~ \J A m\ and m 4 ~ a/A m 2 SND . In this case, we 

can decompose the effective mass in a term well-known from three-flavor analyzes and a 
contribution from the LSND scale m , namely: 


(m) 


3+1B 


(m) 


NH 

3 


+ V^LSND ^1^ 


e41 


where 


(m) 3 = mi cos 2 0 q + a/A 


m 2 & e ia sin 2 


9 e + y Ai/n sill" #chooz • 


(26) 


The term |(m)^ H | corresponds to the effective mass in case of three neutrinos with a normal 


hierarchy. The maximal value of |(m) 3 | in the range of m 2 -C Am 2 is known to be less 
than 0.007 eV [EHE2j, whereas a/A/77 2 snd \U ei \ 2 ~ 0.008 eV and 0.005 eV for the cases 
with A-m 2 g ND = 0.9 eV 2 and 1.8 eV 2 , respectively. Therefore, if mi is negligible and if 
Am 2 SN D = 0-9 gV 2 , there can be no complete cancellation, and a lower (upper) limit on 
(m) of 0.001 (0.015) eV can be expected. If mi is negligible and Am 2 SND = 1.8 eV 2 , 
however, there can be complete cancellation, and the upper limit of (m) is 0.012 eV. For 
both values of Am 2 SND it turns out that also larger values for mi can lead to a vanishing 
effective mass, defining a “cancellation regime” in Fig. 0 If we vary Am 2 SND in the range 
(0.2 - 0.5) eV 2 and mi is small, then (m) is controlled by the term | U e4 1 2 a/ Am 2 snd , which 
is between 0.02 and 0.07 eV, and consequently there is no complete cancellation. For higher 
values of mi there can be complete cancellation in this case also. If the LSND scale is fixed 
to 0.9 or 1.8 eV 2 and mi ^ 0.1 eV, then the masses mi i2i3 are quasi-degenerate and their 
contribution to the effective mass dominates the contribution from 777,4. Consequently, the 
effective mass in this case is (m) 3+1B ~ mo (1 — sin 2 20 0 sin 2 a/2), where mo denotes the 
common mass scale of the three lightest neutrinos and 0 chooz has been neglected. For 
the MiniBooNE range of (0.2 - 0.5) eV 2 this happens for slightly larger values of mi, but 
the important aspect that the effective mass is non-zero in this case holds as well. The 
reason for this is that solar neutrino mixing is non-maximal. All the discussed features are 
reflected in Figure O 

The kinematic neutrino mass is for mi = 0 


m 


3+1B 


~ \! A/77,0 sill' 


00 + sill 2 0 CHOOZ A m\ + |C e4 | 2 Am 2 


LSND 5 


(27) 


which is essentially determined by the term a/A m^ND \U e i\- For Am 2 SND = 0.9 eV 2 as 
well as Am 2 SND = 1.8 eV 2 this product is ~ 0.085 eV, whereas for Am 2 SND in the range 
(0.2 - 0.5) eV 2 it can vary between (0.1 - 0.2) eV. This is reproduced in Figure El It is 
to be noted that m^ +1B can be one order of magnitude larger than the maximal effective 
mass. We can summarize scenario 3+1B for small m 1 as 

£ 3+1B ^ i/Am£ SND > m 3+1B » (m) 3+1B . ( 28 ) 

Looking at Figs. El and El we can make the following statements: if cosmology improves 
the limit on E to be below a J Ami SND , scenario 3+1B can be ruled out. An effective mass 
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Figure 6 : Same as previous Figure for scenario 3+1C. 


above 0.05 (0.01) eV rules out both overlap points at AmL SND = 0-9 and 1-8 eV 2 when 
mi ^ 0.1 (0.01) eV is assumed to be small. In this case a successful KATRIN search will 
make this also possible. 


3.3 Neutrino Masses and Neutrinoless Double Beta Decay in 
Scenario 3+1C 

The mass spectrum for the scenario 3+1C is depicted in Fig. [I] We identify Am|, = m\—ml 
and Am\ = — m\. The heaviest neutrino m 4 is separated by the LSND scale from 

the remaining three, which enjoy an inverted mass ordering. With this identification the 
different masses can be expressed in terms of the smallest mass m\ and the mass squared 
differences as 

m ,2 = \/m\ — A?71q + A m\ , 

m 3 = y/m\ + A m\ , (29) 

m 4 = yj mf + Am 2 SN d • 

In Fig. 01 we show the four masses as well as their sum as a function of the smallest mass 
mi. We always have quasi-degeneracy between m 2 and m 3 , and m\ is quasi-degenerate 
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with them once it is in the vicinity of 0.1 eV. The results on the total sum of masses are 
hardly distinguishable from 3+1B, 


E 3+ c ~ \J Am 2 SND • (30) 

For the mixing matrix elements one finds \U e \\ 2 ~ sin 2 ^chooz as well as 4 \ U ei \ 2 \U^\ 2 = 
sin 2 26 l L g N D and |C/ e o| 2 ~ cos 2 9 e and \U e5 \ 2 ~ sin 2 # 0 . 

We show in Fig. El our result for the effective mass as well as for the kinematic mass as a 
function of the smallest mass. The effective mass in the limit of small m\ is 


(m) 


3+1C 


a/ Am 


0 ij3 


COS' 


e Q 


+ e ia sin 2 6, 


o 


) + 


I U ( 


e4 


(31) 


which can also be written as 

v 3+1C 


/ \ C 

(m) 


(m)' H + J A m£ SND e 17 \U e ,\ 


(32) 


where the term Km)™) corresponds to the effective mass in case of three neutrinos with 
an inverted hierarchy [33, • The contribution of the smallest mass m 4 plays a sub-leading 
role as it is also multiplied by small sin 2 #chooz- The absolute value of (m)™ is known 
to be y/Am\y/l — sin 2 29 & sin 2 (/3 — a )/2 and lies (for m\ -C A m 2 A ) between 0.007 and 
0.05 eV <Tl. i(T?l . Whether there is complete cancellation in (m) 3+1C or not depends on the 
value of a/A raLSND ItT e 4 1 J • For Am 2 SND = 0.9 (1.8) eV 2 it is given by ~ 0.008 (0.005) eV. 
Therefore, for AmL SND = 0.9 eV 2 there can be complete cancellation as is seen in the first 
panel of Fig. El but since for Am 2 SND = 1.8 eV 2 the value of a/A/Ulsnd £4 4 1 2 is less than 
0.007 eV, there is no complete cancellation in this case (note that it is the other way around 
in scenario 3+1B). The lower limit on the effective mass is roughly 10 -3 eV. Similarly, in 
the case where we vary A'/n 2 gND from (0.2 - 0.5) eV 2 and |t/ e 4 1 2 from 0.05 to 0.1, the value 
of the product a/A m 2 SND |f / e4 | 2 can be smaller than 0.06 eV and therefore in this case also 
there can be complete cancellation over a wide range of rri\. The maximum value of the 
effective mass in the small mi limit is roughly 

- V Am A + V^LSND WeA? ■ (33) 

Since the product a/A m 2 SND |C/ e 4 1 2 is approximately the same for 0.9 and 1.8 eV 2 , it follows 
that ( m fakes the same value (~ 0.07 eV) in these two cases. If rn i is around 0.1 eV, 
then the three lightest masses are quasi-degenerate, and similar comments as for scenario 
3+1B discussed in the previous Subsection apply. 

The kinematic neutrino mass for m 4 ~ 0 is 

mj +ic ~ ^A m\ + \ U e4 \ 2 Aml SND , (34) 

which is somewhat larger than the maximum effective mass. The blue (dark) band in 
Figure El shows mp against m \. ft is to be noted that in the first two columns of Figure El 
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the width of mp is due to the variation over the allowed range of A m 2 A as Awl SND is held 
fixed in these plots. In the third column the width is due to variation of both A?77 2 SND 
and A m\. Let us summarize the situation: 

£ 3+1C - \j Am£ S ND > m l +1C ~ • ( 35 ) 

Ruling out scenario 3+1C could be achieved if cosmology improves the limit on E below 
a JAm‘l s - ND . For small rri\ < 0.1 eV and a successful KATRIN search both overlap points 
at At? 1 l SND = 0.9 and 1.8 eV 2 are ruled out. An effective mass above 0.07 eV rules out all 
three cases under discussion, unless mi ^ 0.1 eV. 


4 Four Neutrino Mass Matrices 

Models incorporating an extra sterile neutrino have been developed in many papers m 
EH ESI ESI El!- In this Section we wish to summarize the typical mass matrices that are 
consistent with the experimental data in 3+1 scenarios. In particular, we look for simple 
U( 1) flavor symmetries which can force the approximate form of the mass matrices. The 
mass matrices for the 2+2 scenarios are discussed in Appendix iBi 

3+1 scenarios have the property that the sterile neutrino does practically not participate in 
solar and atmospheric neutrino oscillations. Consequently, there is very little dependence 
on the respective sterile neutrino fraction. Another general aspect of the results in 3+1 
scenarios is that the well-known three-flavor mass and mixing matrices (see the overviews 
hi jjSBI E2]) are “embedded” in the four-flavor mass and mixing matrices. This means, 
in particular, that the corrections to the usual three-flavor mass matrix are of order 5 
?7i4 sin 2 6 l LSND ~ A 2 a/Aiblsnd — We introduced here a small parameter A ~ 0.1, 

to estimate the different mass and mixing scales in the four neutrino framework. Both the 
LSND and the CHOOZ mixing angle are assumed to be of order A, and the mass scales 
are related through A m 2 & — A 2 A m\ ~ A 4 A/Blsnd- I n the approximation we are using, 
terms of order A 2 a/Ati7l SND are subleading. 

4.1 The Mass Matrix in Scenarios 3+lAa and 3+lAb 

In case of scheme 3+lAa, one has 

A cos 9 & sin # 0 \ 

sin 6a — sin#© cos 6a cos#© cos 9a p 
cos 9a sin 9 & sin 9a — cos # Q sin 9a 1 ' 

A A A / 

where we set the Dirac phases to zero and included small terms of order A ~ 0.1 without 
writing possible order one coefficients. These terms indicate the typical order of both the 

5 We have been made aware of an upcoming analysis on this subject EE- 


3+lAa ^ 


( A 
A 
A 
1 
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CHOOZ angle sin Whooz and the LSND parameter sin $lsnd• The above mixing matrix 
is unitary only to order A. In principle, the order one entries receive additional terms of 
order A to cure this. The following analysis, however, is not harmed by this. With a given 
mass hierarchy we can obtain now the approximate form of the mass matrix. By looking 
at Fig. 01 we can see that typically m 4 ~ m 3 ~ m 2 mi holds. So, setting m\ = 0, we 
have 

c A c e s Q (e l/3 - e* 7 ) SAC 0 s 0 (e l/3 - e* 7 ) A ^ 

e ia s\ + c\(e l P Sq + e 11 cfy CAS\(e ia — e^CQ — e l @ Sq) A 

e ia c\ + s\(e z P Sq + e 1 ' 7 c@) A 

0 

Here we have defined ca = cos 9 a and sa = sin 9 a- A matrix with the entries of the s 
column zero and the remaining elements of order one conserves the flavor charge L s . 
Obviously, the upper left 3 x 3-block of m v corresponds to the well-known three flavor 
mass matrix in case of quasi-degenerate neutrinos. Apart from the usual /i-r exchange 
symmetry Ena. we can have several interesting special cases: depending on the relative CP 
parities of the three heavy neutrinos, and setting for simplicity 9 a = 7 t /4 and 6 q = tt/ 4, 
we can have 3x3 matrices proportional to the unit matrix (a = (3 = 7 = 0) or with only 
a non-vanishing ee and fir element {(3 = 7 = tt and a = 0 ) fUBJ. 

The mixing matrix in scenario 3+lAb is obtained by exchanging the second and fourth 
row of the mixing matrix, i.e., 




m 3+ 1A a _ J A m 2g ND 


f/3+lAb ^ 


( A 

A 

A 


V 1 


sin 9 q 

cos 9 q cos 9 a 
— cos 9 q sin 9 a 
A 


cos 9 q A \ 

— sin 9 q cos 9 a sin 9 a 
sin# 0 sin^A cos 9 a 
A A / 


(37) 


The mass matrix looks identical to case 3+lAa, the only change being the replacement 
a 7 . This is analogous to the three-flavor case, in which the structure of the mass matrix 
for quasi-degenerate neutrinos does not depend on whether the neutrinos are normally and 
inversely ordered. 


4.2 The Mass Matrix in Scenario 3+1B 

The mixing matrix is given by 


^y3+lB 


/ cos 9 q 
— sin 9 q cos 9 a 
sin0 0 sin#4 

V A 


sin 9 q A 

cos# 0 cos 9 a sin^A 
— cos 9 q sin 9 a cos 9 a 
A A 


A \ 

A 

A 

1 / 


(38) 


Regarding the mass states, we have for a smallest mass m\ ^ 0.005 eV that m 4 ~ 
a/A?71lsnd ^ m 3 — a/A m\ ~ m 4 A m 2 — a/A m\ ~ m 4 A 2 . The mass matrix 


18 












then reads 


m 3+1B ~ 


^ m LSND 


/ 0 0 0 

• s \A C A S A A 

• • 4A 

V • • 


A 

A 

A 

1 


\ 


Again, the upper left 3 x 3-block corresponds to the well-known three-flavor mass matrix, 
which can be obtained by demanding L e to be conserved 6 . A four-flavor mass matrix with 
only the ss element non-zero conserves 7 L e + + L T . As long as rn\ (and therefore also 

m 2 and m 3 ) are one order of magnitude below 7714 , the main structure of the mass matrix 
remains, i.e., the ss entry is one order of magnitude larger than the remaining ones. 


4.3 The Mass Matrix in Scenario 3+1C 

The mixing matrix is given by 


f/3+lC 


( A 

sin 9 a 

cos 9 a 

V A 


cos 0 © 

— sin 0 © cos 9 a 
sin sin 6a 
A 


sin 0 q 

cos 0 Q cos 9a 
— cos 0 Q sin 9a 
A 


A \ 

A 

A 

1 / 


(39) 


The three light neutrinos correspond approximately to the well-known inverted hierarchy 
case of three neutrinos. With mi ~ 0 and m 3 ~ m 2 — \f A m\ ~ m 4 A ~ \J Am^ A, we 
get 


m 


3+1C 

V 





2 

LSND 


( e ia CQ + Sq A ( e ia — e*0)cAS©c© A 
(e* a s© + e*0 Cq)c\ A 


( e ia _ e *0)sAS©c o A A ^ 
(e m s| + e*0 c|)c A -sa A A 
(e ia SQ +e 10 Cq)s 2 a \ A 


V 



The full four-flavor matrix conserves approximately (i.e., when we neglect A) the flavor 
charge L e + + L r . We have again for the upper left 3x3 block the well-known three 

flavor mass matrix of an inverted hierarchy, which displays for 0a = tt/ 4 a /i-r symmetry. 
If we set 0© = 7 t/4 and choose 8 a = 0 and (3 — it, then all entries except the ss, e/i and er 

6 Conservation of L e in case of four flavors would be achieved when the /is and rs entries are also of 
order one. 

7 Obviously, conserving L e + L M + L T in case of three flavors leads to Dirac neutrinos. For four flavor 
scenarios one would have to ask for L e + L M + L T + L s conservation. 

8 Note that such a Pseudo-Dirac structure will lead to enhanced stability with respect to radiative 
corrections. 
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elements vanish to order A: 


m 3+1G ~ 


A' m LSND 


/ 0 1 1 0 \ 

• 0 0 0 

•■00 


V 


1 / 


(40) 


The global symmetry forcing this form of the mass matrix is L e — — L T , which was 

introduced hrst for the three-flavor case EU, but was used also for the four flavor case fST]. 


5 Comments on 3+2 Scenarios 

By means of introducing more sterile neutrinos, the goodness of fit for explaining the 
LSND and other short-baseline data can, not really surprisingly, be improved. In this 
respect, the 3+2 scenario has been put forward to make the interpretation of all neutrino 
data less problematic m- These schemes have three neutrinos actively oscillating among 
themselves and two additional sterile neutrinos responsible for the LSND anomaly. Models 
to accommodate 3+2 scenarios can be found in :72j- With 5 neutrinos participating in 
neutrino oscillations, 4 independent Am 2 are present. In addition to the three discussed 
previously, we have to deal in addition with Am 2 ]. 

In the analysis of Ref. HE! two best-fit points are given, one of which corresponds to 
Am| x = 22 eV 2 , which is clearly not consistent with cosmological constraints. The second 
best-fit value is at 

A = 0.46 eV 2 , A ml, = 0.89 eV 2 , 

,51 ( 41 ) 

U eA = 0.090 , U e5 = 0.125 , 4 = 0.226 , U^ = 0.160 . 

This identification of the mixing matrix elements assumes that the three active neutrinos 
are lighter than the two sterile ones, i.e., a situation resembling scenarios 3+1B and 3+1C. 
We can decompose the effective mass as a term from the three active neutrinos and a term 
from the two sterile ones, i.e., 

(m) 3+2 = |(m) 3ac + (m) 2st | . (42) 

In this case, the two additional mass scales imply an additional contribution to the effective 
mass, reading 

(m) 2st = \U eA \ 2 m 4 e* 7 + |t/ e5 | 2 m 5 e* p = \U e4 \ 2 \Jm\ + Amh e* 7 + \U e5 \ 2 \Jm\ + Am|, e ip ,(43) 

where we have introduced a fourth relevant Majorana phase p. If mi = 0, then the best-fit 
values given above yield |(m) 2st | ~ (0.01 - 0.02) eV, where the range is caused by the effect 
of the Majorana phases. 
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The cosmological mass parameter is E ~ \Jm\ + Am^ + m\ + A m 2 ± , which for mi = 0 
is 1.6 eV, remarkably close to the current relevant limit from cosmology (1.64 eV) obtained 
in Ref. |2B| - 

The contribution to the neutrino mass measurable in KATRIN is roughly 

mg ~ \Z\U e4 \ 2 (ml + A m%) + \U e5 \ 2 (mj + Am|) (44) 

which for mi = 0 is 0.1 eV, but can reach testable values if mi ^0.1 eV. 

If the three active neutrinos display an inverted hierarchy, then their contribution |(m) 3ac | 
to the effective mass (see scenario 3+1C in Section lies between (m) = \JAm\ and 
(m) = > J A m\ cos 20 q, or numerically: (0.007-0.06) eV. Complete cancellation between 
(m) 3ac and (m) 2ht is possible. On the other hand, positive interference of the two terms 
can lead to (m) ~ 0.08 eV. If the three active neutrinos enjoy a normal mass ordering 
with mi = 0, then their contribution (m) 3ac to (m) does not exceed 0.007 eV (see scenario 
3+1B in Section 13.21) , so that (m) 2ht dominates and represents basically the prediction for 
(m). Note however that these considerations take use of the best-fit values of the fourth 
and fifth neutrino sector. Varying these parameters within their allowed range and/or 
effects of non-zero m\ might easily allow for severe cancellation. Moreover, one might also 
exchange Am| : and Am^, leading to identical values for E, slightly larger values of mg 
and |(m) 2hi | ~ (0.003 — 0.02) eV (all if mi = 0). Consequently, the effective mass vanishes 
for both orderings of the three active neutrinos. The upper limit is 0.02 (0.08) eV for 
normally (inversely) ordered active neutrinos. 

Up to now the two sterile neutrinos were assumed be heavier than the active ones. Also 
possible is that the two additional neutrinos are lighter than the three active ones, thereby 
resembling scenarios 3+lAa and 3+lAb. For the mixing matrix elements one has to ex¬ 
change the indices 4 <->■ 1 and 5 2. The three active ones generate an effective mass 

larger than roughly a/ Am^ + m\ cos 26 Q , which for mi = 0 is 0.1 eV (see Section mi) . 
The contribution of the additional neutrinos is significantly suppressed. The KATRIN pa¬ 
rameter is approximately a/A m^ ~ 0.9 eV, surely a testable value. Cosmology will have 
to probe E ~ 3 a/A m 2 -^ ~ 2.8 eV, which is incompatible already with current limits. 

Finally, it should be clear that the known three-flavor mass and mixing matrices are em¬ 
bedded in the 5x5 mass and mixing matrices. We note that scenarios with 3 or more 
sterile neutrinos will also show this “embedding”. 

6 Conclusions and Summary 

We have examined the constraints on LSND induced scenarios with extra sterile neutrinos 
from current and future bounds on neutrinoless double beta decay, tritium beta decays 
and cosmological limits on the sum of neutrino masses. Since 2+2 scenarios are already 
disfavored by the present solar and atmospheric data we considered the 3+1 scenario in 
the main part of the paper. The values of Am 2 SND considered in our analysis are 0.9 eV 2 
and 1.8 eV 2 , allowed by a combined analysis of SBL + atmospheric and K2K data mm3 
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We also considered Aml SND to vary in the range (0.2 - 0.5) eV 2 , motivated by MiniBooNE 
sensitivity plots. For sake of completeness we discuss the 2+2 scenarios in Appendix iBl 
Within the 3+1 scenario there are three possibilities, 3+lAa,b, 3+1B and 3+1C. The sum 
of neutrino masses is ~ 3(1, l)y/A?n 2 SND , for scenario 3+lAa,b (3+1B, 3+1C) neglecting 
the solar and atmospheric mass differences. The effective mass can be written as a known 
three-flavor contribution plus an additional term stemming from the LSND scale. Let us 
summarize the different aspects: 

• in the 3+lAa and 3+lAb scenarios we have three quasi-degenerate neutrinos of mass 
a/ AmL SND and a fourth state separated from them by the LSND mass scale. The 
sum of the masses is ~ 3A/Am 2 SND in the limit where the mass of the lowest state 
is vanishingly small. This implies small values of A/n 2 SND (~ 0.1 eV 2 ) to comply 
with the cosmological bound of order 1 eV for the sum total of the masses. For the 
considered values of Am 2 SND the sum of masses already exceeds the cosmological 
limit of ~ 1 eV. Therefore this scenario is highly constrained from cosmology in 
particular for higher values of AmL SND . The effective mass measured in neutrinoless 
double beta decay varies from a/Aj?Il SND cos 2 9 & to Am 2 SND . Since 6 & = 7 t/4 is 
no longer allowed by the current data, the minimum value of the effective mass is 
non-zero in this case 9 . Depending on the value of Am 2 SND , the lower value of (m) 
varies from ~ (0.04 — 0 . 2 ) eV. A part of this region is already disfavored by the limit 
from the Heidelberg-Moscow experiment. Low values of sin 2 9q and large values of 
Am 2 S ND jeopardize scenarios 3+lAa and 3+lAb because the minimal (m) becomes 
too large. The neutrino mass measured in beta decay in this scenario is ~ \J AmL SND 
and coincides with the upper limit of the effective mass in neutrinoless double beta 
decay; 

• the 3+1B scenario corresponds to three neutrino states with a normal hierarchy 
separated from a fourth state by the LSND mass scale. For small mi < 0.01 eV the 
sum of the masses is given by S 3+1B ~ a/A m 2 SND . Hence, the cases A'/Ulsnd = 0-9 
eV 2 and Am 2 SND in the range (0.2 - 0.5) eV 2 are consistent with the cosmological 
mass bound. Whether there is complete cancellation for the effective mass depends 
upon the product |D e 4 | 2 A/Am 2 SND . For our choice of parameters and hierarchical 
masses we got very low values of (m) only if Am 2 SND = 1.8 eV 2 . The mass measured 
in tritium beta decay is larger than (m) for rri\ ^0.1 eV; 

• scenario 3+1C corresponds to the usual three generation inverted hierarchy picture 
plus an additional neutrino at a higher scale separated by the LSND gap. The 
second and the third state are quasi-degenerate at \/ A m\. For mi > 0.1 eV they 
become quasi-degenerate with the first state. The cosmological constraint on the 
3+1C scenario is approximately the same as for 3+1B. The prediction for (m) can 
be split, into a contribution from a three neutrino inverted hierarchy scenario plus 
a term e* 7 U e 4 | 2 \J Am 2 SND . The contribution of the three-flavor term to (m) is 

9 This also is the case in scenario 2+2A. 
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between (0.007 - 0.05) eV. Hence, in this case we get very low values of (m) for 
A?71 lsnd = 0-9 eV 2 and in the range (0.2 - 0.5) eV 2 . 

Regarding the four neutrino mass and mixing matrices, 3+1 scenarios “embed” in their 
mass and mixing matrices the well-known three-flavor matrices. Corrections to these three- 
flavor matrices are of order a/A/Ulsnd suA^lsnd, he., of order Am 0 , for the mass matrix 
and of order sin $lsnd for the mixing matrix. The 3+1B and 3+1C scenarios can in 
principle be motivated by an approximate L e +L^+L r global symmetry 10 , whereas scenarios 
3+lAa,b correspond to a L s global symmetry. Scenario 3+1C can also be motivated by 
conservation of L e — — L r . 
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Note added : The recent results from three year WMAP data in combination with large- 
scale structure and supernova data give a bound on the sum of three neutrino masses 
to be Hm v < 0.68 eV (95% C.L.) which is not very different from the first year limit 
EH- Therefore the bound on the sum of neutrino masses that we use in this paper is 
not expected to change significantly with this new data. More stringent bounds can be 
obtained by including the Lyman-ct forest data m- 


A Oscillation Probabilities for Short Baseline Expe¬ 
riments 


Here we list the relevant probabilities for short-baseline experiments in the 2+2 and 3+1 
scenarios. The most general expression for neutrino survival or conversion probability for 
N neutrino generations is given by 


Pv a vp — fia/3 


4 ^2 Pod U/3i U aj U Pj sin 2 j , 

j>i ' ' 


(45) 


where i,j varies from 1 to A for N generations, and \j = 2.47 (E u / MeV) (eV 2 /Am 2 ) 
m. If the Ami corresponds to A m^, A m\ or Am 2 SND , we denote it with A 0 , Aa and 
-^lsnd; respectively. The actual form of the various survival and transition probabilities 

10 The 2+2 scenarios, discussed in the Appendix, have other candidates for global symmetries, for instance 
L e + — L t — L s . Moreover, there can be a r-s exchange symmetry in analogy to the successful /i-t 

exchange symmetry for three neutrinos. 
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will depend on the spectrum of Am 2 chosen. The above Eq. © is assuming the Ce¬ 
phases to be zero. Below we list the relevant probabilities for short-baseline experiments 
in the 2+2 and 3+1 scenarios for the reactor experiments Bugey [Q|, CHOOZ 0 and the 
accelerator experiments CDHSW ljj . LSND [1] and MiniBooNE [53] . The energy and 
length scales involved are such that one mass scale dominance approximation holds true 
in most of the following cases (excepting CHOOZ): 

• Scenario 3+lAa: 

Pv e v e = 1-4 O 2 (1 - U 2 el ) sin 2 (7rL/A LSND ) (Bugey) 

P VeVe = 1 - 2 U 2 el (1 - U 2 el ) - 4 U 2 e2 sin 2 (7rL/A A ) (CHOOZ) 

JVm = 1 - 4 U h (! - u h) sin 2 (vrL/A LS ND) (CDHSW) 

Pv,„v K = 4 U 2 el U 2 A sin 2 (7rL/A LSN D) (LSND, MiniBooNE) 

The relevant formulae for scenario 3+lAb are obtained from the above by replacing 
2 4. 

• Scenario 3+1B: 


Pv.v, = 1 - 4 lift (1 - (A. 2 ,) sin 2 (iri/A LSND ) (Bugey) 

Pv.v, = 1-2 t/ 2 „ (1 - t/ 2 „) - 4 t/ 2 3 siu 2 (yL/A A ) (CHOOZ) (46) 

P^v r = 1-4 1% (1 - 1%) sin^L/Assi) (CDHSW) 

= 4 C/ 2 4 1% sin 2 (jrL/A LSND ) (LSND. MiniBooNE) 

Note that the probabilities in the 3+1B picture can be obtained from 3+lAa by replacing 
U e i by U e 4, by and U e2 by U e3 in the CHOOZ probability. The probabilities for 

the scenario 3+1C are the same as those of 3+1B excepting U 2 3 in CHOOZ probability 
is to be replaced by U 2 A . In the presence of CP violation the mixing matrix elements are 
complex. In the above 3+1 oscillation probabilities one then has to replace U 2 { with \U a i\ 2 . 

Let us for the sake of completeness also give the oscillation probabilities in the 2+2 schemes: 
• Scenario 2+2A: 

Pv.v, = 1-4 (C/ 2 , + Ul 2 ) (1 - Ul - f4) sin 2 (7rL/A LSND ) (Bugey) 
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Pp.v. = 1 - 2 (E& + f/ 2 2 ) (1 - U 2 el - U 2 e2 ) - 4 C& P 2 2 sin 2 (7rL/A A ) (CHOOZ) 


P VeVe = 1-4 (C/ 2 , + P 2 2 ) (1 - C/ 2 , - f/ 2 2 ) sin 2 (vrL/A L sND) (CDHSW) 

Pv,ve = 4 (C/ e i i + P e2 C/ m2 ) 2 sin 2 (7rL/A LS ND) (LSND, MiniBooNE) 

The survival and oscillation probabilities in scenario 2+2B are obtained from those of 
scenario 2+2A by making the change 1 3 and 2 <->• 4. In case there are Dirac phases 

in the PMNS matrix, again U 2 % will have to be replaced with \U Q i\ 2 . In what regards the 
(LSND, MiniBooNE) probability, the relevant term would read U* A U fl \ + U* 2 U ll2 instead 
of (Pei Up^UaUpz). 
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Figure 7: All four neutrino masses and their sum E against the smallest neutrino mass m\ 
for scenario 2+2A (left) and 2+2B (right). The masses m3 and m 4 are not distinguishable 
in the plot. We choose Am 2 SND = 0.2 eV 2 . 


B Neutrino Masses and Neutrinoless Double Beta 
Decay in 2+2 Scenarios 

The 2+2 scenarios are disfavored by the combination of solar and atmospheric data re¬ 
gardless if LSND results are confirmed by MiniBooNE or not. However, for the sake of 
completeness we discuss in this Section the constraints on 2+2 mass spectra from cosmol¬ 
ogy, neutrinoless double beta decay and tritium beta decay. 


B.l Neutrino Masses and Neutrinoless Double Beta Decay in 
Scenario 2+2A 


On the left side of Fig. 0 the mass ordering of scheme 2+2A can be seen. In this scheme 
it holds that Am|, = Am 43 and A m\ = Am^. One has 

m 2 = a Jm\ + A m\ , 

m 3 = yjm\ + Am 2 SND - Ara| , (47) 

m 4 = sjm\ + Am 2 SND . 

The left side of Fig. 0 shows the four mass values and their sum E as a function of the 
smallest mass m Due to AmL SND A m\ the masses m 4 and m 3 are hardly distin¬ 
guishable in the plot and for rri i ^ 0.1 eV the two lightest masses rri\ and m 2 become 
quasi-degenerate. To clarify the role of the cosmology bounds on E, we included in the 
Figure a value of E = 1 eV. To let the sum of neutrino masses lie below this limit, small 
values of Am 2 SND and of rn i are implied: for small mi ^ 0.05 eV it holds that 


m 4 ~ ?« 3 ~ 



LSND ^ m 2,l 


(48) 
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Figure 8: The effective mass (m) and the mass mp measured in beta decay experiments 
plotted as a function of the lowest mass in 2+2A scenario. The left hand column is for 
AmL S ND = 0.3 eV 2 and U 2 { = f/ 2 2 = 0.002. The right hand column is for Am-LSND = 0.9 
eV 2 and U 2 X = f/ 2 2 = 0.005. The other parameters are varied in their current 3cr allowed 
range and all the phases are varied between 0 and 27 t. Also shown is the mass mp that 
will be measured in beta decay experiments, the current and a prospective future limit on 
the effective mass and the future KATRIN limit. 


Hence, £ ~ 2^AmLg ND and therefore low values of A?71 l SND ^0.3 eV 2 are implied by the 
condition £ ^ 1 eV. 

Turning to the constraints on the mixing matrix elements, \U e \\ 2 and |t/ e 2 1 2 are constrained 
by the short baseline reactor experiment Bugey and the reactor experiment CHOOZ. 
sin 2 20 L g N D is given by the combination {U eA U^\ + U&U^) 2 and the combination {\U^i\ 2 + 
| U^ 2 1 2 ) will be constrained by the CDHS experiment. We give the relevant expressions for 
the probability in Appendix El Since the one mass scale dominance approximation holds, 
one can use two-parameter plots to find the constraints on these mixing parameters. How¬ 
ever, for the sake of illustration in this paper we use the values of AmLg ND and sin 2 2$lsnd 
from the MiniBooNE sensitivity plot given, e.g., in 0. We take the representative val- 
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ues 11 for (A/?t 2 snd , sin 2 20LSND) as (0.3,0.02) and (0.9,0.008). To extract U e \ and U e 2 from 
sin 2 2#lsnd we make the plausible assumption U e \ — U e 2 and U ’ 2 X = U 2 2 = 0.5 as implied 
by atmospheric data. This assumption was for instance used in [,4Ej • With this assumption 
we have f/ 2 , = U 2 2 = 0.002 (0.005) for Am 2 SNn = 0.9 (0.3) eV 2 . In any case it is to be 
noted that in the 2 + 2 A scenario U e 1 and U e 2 multiply the smaller masses uq and m 2 , and 
as we will see below their contribution to effective mass as well as the mass measured in 
beta decay is sub-leading. We furthermore have \U e s\ — cos0 o and | U e 4 ~ sin 0 Q . 

Since mi and m 2 are small in scenario 2 + 2 A and in addition multiplied with the small 
elements \U e i\ 2 and \U e2 \ 2 , respectively, we can neglect terms including these quantities in 
what follows. Then the effective mass in scenario 2+2A reads 

(m) 2+2A ~ yj A ml SND yjl- sin 2 20 0 sin 2 (/ 3 - 7)/2 . (49) 

The non-maximality of solar neutrino mixing implies therefore a non-vanishing effective 
mass. This is in analogy to the three-flavor case with an inverted hierarchy or quasi¬ 
degenerate neutrinos. Choosing for instance \J ~ 0.5 eV and the values of 0 0 
from Eq. we can predict that 

cos 200 yj AmL SND ~ 0.1 eV k) (m) 2+2A k) yj Am 2 SND k) 0.5 eV . (50) 

This range of (m) is well within reach of currently running or planned 0 u/3/3 experiments. 
In Fig.|Hlwe show the effective mass as a function of the smallest neutrino mass mi. We also 
show the present bound from the Heidelberg-Moscow experiment in this Figure, together 
with a prospective future limit. It is to be noted that some part of the regions are already 
disfavored by the Heidelberg-Moscow limit. Thus non-maximality of solar neutrino mixing 
angle coupled with the existing limit from Heidelberg-Moscow experiment already puts 
some constraint on the 2 + 2 A mass pattern. 

Neglecting terms proportional to U^mi and C 2 2 m 2 , we have for the kinematic neutrino 
mass 


mJ +2A ~ y/\U e3 \ 2 ml + \U ei \ 2 m\ ~ yj Am£ SND . (51) 

Since cosmology implies that \J Am 2 SND is below roughly 0.3 eV, we expect that m^ should 
be close to the lowest value reachable by KATRIN, but close to the current limit on the 
sum of neutrino masses from cosmology. Since 

m l +2A ~ \j Am LSND - E 2+2A /2 ~ yj Amj( SND yjl — sin 2 20 o sin 2 (/5-7)/2 , (52) 

one can in principle obtain a set of consistency checks of scenario 2 + 2 A, and obtain some 
information on the Majorana phase combination (3 — 7 . 

n It is to be noted that the constraints on mixing angles from SBL experiments are not as severe in the 
2+2 case as in the 3+1 case. 
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B.2 Neutrino Masses and Neutrinoless Double Beta Decay in 
Scenario 2+2B 

The right panel of Fig. [3 shows the values of the four neutrino masses and their sum S in 
scenario 2+2B, in which one has A m 2 Q = m\ — ml an d = m \ ~ m 3 - We can write the 
masses of the neutrino states in terms of the smallest mass m 4 and the three mass squared 
differences as 

m 2 = \Jml + A m 2 Q , 

m 3 = y/ml + Am£ SND - A m\ , (53) 

m 4 = \Jml + Am LSND • 

For the mass values similar statements as for scenario 2+2A apply, namely that m 4 ~ m 3 ~ 
a/A mL SND ~ S 2+2B /2, where Am^ND should lie at the lower end of its allowed range in 
order to obey the constraints from cosmology. The smallest neutrino mass m 4 should lie 
below 0.05 eV and therefore m 2 and m 4 are typically not very close to each other, unless 
mi ~ 0.02 — 0.05 eV. 

In what regards the mixing matrix elements, scenario 2+2B is obtained from scenario 
2+2A by exchanging the indices 1 3 and 2 4. Hence, \U eA \ and \U e2 \ are roughly 

given by cos^q and sind Q , respectively. The elements \U e3 \ and U e4 \, however, are implied 
to be small. 

We show in Fig. |21 the effective mass as a function of the smallest mass. The effective mass 
is approximately given by 


(m) 


2+2B ~ |cos 2 6q m\ + sin 2 d Q m 2 e ia + \U e3 \~ m 3 e 10 + \U e4 \ 2 m 4 e* 7 
cos 2 9 e mi + sin 2 9 e m 2 e ia + a/A m 2 snd (|H e3 | 2 e 10 + \U eA \ 2 e* 7 ) 
sin 2 d Q y/Am 2 Q + y/ A ml SND (| U e3 \ 2 e^ 0 - a) + | U e4 \ 2 e i(7_/3) ) , 


(54) 


where we neglected mi for the last approximation. 

Since U 2 3 and U 2 A are small, the two large masses m 3 and m 4 are multiplied with small 
mixing matrix elements, whereas the small masses ni\ and m 2 are multiplied with large 
mixing matrix elements. As a consequence, there can be cancellations leading to a very 
small or zero effective mass. Note the analogy of this situation with the three-flavor case: in 
the inverted hierarchy the large masses are multiplied with mixing matrix elements corre¬ 
sponding to the large solar neutrino mixing, whose non-maximality allows no cancellation. 
In the normal hierarchy, the largest mass m 3 is multiplied with the smallest mixing matrix 
element, and complete cancellation can occur. Since the degree of cancellation depends on 
the values of the two small quantities U 2 3 and I/ 2 4 for this case, we present our results for 
Aut-lsno = 0-3 eV 2 and two different choices for U 2 3 and I/ 2 4 . In the left panel we show 
the plot where we assume U 2 3 ~ H 2 4 = 0.005. In this case complete cancellation can occur. 
In the right panel we present our results with U 2 4 = 0 and U 2 3 = 0.01. In this case the 
minimal (m) has a higher value. The two terms in the last equation have the typical order 
a/ Am 2 sin 2 9 e ~ (2 — 4) • 10 -3 eV and a/ A/ttlsndI^ e31 2 — 0.005 eV. The upper limit on 
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Figure 9: Same as previous Figure for scenario 2+2B and Awlsnei = 0-3 eV 2 . The left 
hand column is for I7 2 3 = U 2 4 = 0.005, the right hand column is for U 2 4 = 0 and U 2 3 = 0.01. 


(m) 2+2B is then roughly given by 0.01 eV. Hence, we can in principle distinguish scenario 
2+2A from 2+2B via 0 u/3/3 as long as mi is small. This is analogous to the situation 
normal vs. inverted hierarchy in the three flavor case. 

The kinematic neutrino mass as measurable in the KATRIN experiment is given by 

m 2 p +2B = ^/\U el \ 2 m 2 i + |R e2 | 2 m 2 + \U e3 \ 2 mj + \ U eA \ 2 m 2 4 
— \/cos 2 6> 0 ml + sin 2 6> 0 m 2 2 + Am£ SND (\U e3 \ 2 + \U e4 \ 2 ) ( 55 ) 

- ^/Am| sin 2 0 0 + A/+ SND (\U e3 \ 2 + \U e4 \ 2 ) , 

where we neglected again m\ in the last approximation. Both terms are of similar magni¬ 
tude and we can expect that ?+ +2B ~ 0.1 eV, larger than the effective mass by an order 
of magnitude and below the future KATRIN limit. 


C The Mass Matrix in the 2+2 Scenarios 

Now we discuss the form of the mass matrices in 2+2 scenarios. Our approach and the 
approximations made are the same as for the 3+1 case, and details are given in Section HJ 
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C.l The Mass Matrix in Scenario 2+2 A 


For scenario 2+2A we can express the mixing matrix as follows EaiBj: 


JJ2+-2A ^ 


( A 

COS 9 atm 

- sin r) sin 6 > atm 
\ - cos Tj sill flatm 


A 

Sin 6> at m 

sin rj COS 6> at m 
COS T\ COS 0 a tm 


cos 9 q 
A 

— cos T) sin 9 & 
sin rj sin 9q 


sin 9q \ 
A 

COST] COS Oq 

— sin ri cos 9 e y 


( 56 ) 


The parameter rj indicates inasmuch sterile neutrinos participate in atmospheric or solar 
neutrino oscillations. For rj — 0 atmospheric neutrinos oscillate completely into sterile ones 
and for rj = 7 t /2 solar neutrinos oscillate into sterile ones. With a given mass hierarchy we 
can obtain now the approximate form of the mass matrix. Glancing at Fig. [3 we identify 
two interesting possibilities, namely 


(i) a/A mL SND ~ m A ~ m 3 » m 2 ~ v^" 1 © ^ m i > 

(ii) a/A mL SND — m 4 ~ m 3 m 2 ~ ~ 0.1 eV . 


(57) 


The first case (i) corresponds to a very small mass rri i and the second one (ii) to two 
quasi-degenerate pairs, though only a small range of mi values allows for this possibility. 
Since the form of m„ is similar in both cases, we mainly discuss case (i). We have then 
m 2 — 772,4 A 2 and the mass matrix reads 


/ e 


if) „2 


772 


2+2A 


^ m LSND 


Cq + e* 7 , 


’© 


V 


A c v (e * 7 — e l/3 ) 

0 A 

• c 2 (e i7 4 + e^ 4 ) 


s v (e* 7 — e l/3 ) \ 

s v A 

(e* 7 4 + e 10 4) 
s 2 (e* 7 4 + e i/3 4 ) / 


We defined the obvious notation c Q = cos # 0 and s Q = sin 0 Q . Terms of order A 2 and 
unimportant factors, such as a coefficient c 4 e* 7 c Q — c 2 e 10 s & for the /it element, are not 
included in our expressions. The factors c lj2 depend on the precise values of the CHOOZ 
and the LSND angles. Note that contributions of the atmospheric mixing are suppressed 
in the mass matrix. If the heavy states m 3 and 7774 have equal CP parities, or when (3 — 7 , 
then this leads to the vanishing of the er and es entries of m u , independent of rj. In case of 
opposite CP parities of z / 3 and u A (which would imply enhanced stability with respect to 
radiative corrections), the ee element and the ts block of m u would be slightly suppressed 
by a factor cos 2 # 0 . Assuming as yet another approximation that 6@ = 7 t /4 would make 
these entries vanish. Several special cases can be obtained from the above matrix. For 
instance, if solar neutrinos oscillate entirely in sterile neutrinos, i.e., rj = 7t/ 2, then 


/ e 10 4 + e 
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777 2+2A ~ 


^ m LSND 


0 0 
0 0 

■ 0 


3*7 


0 if3 


\ 


0 

0 


e i7 4 + e 


if) „2 


(59) 
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This matrix conserves the flavor charge L^ + L T . If we choose equal CP parities of zz 3 and 
z/ 4 , or when fi — 7 , then the er entry vanishes. 

In analogy, if atmospheric neutrinos oscillate entirely in sterile neutrinos (7 = 0) then one 
finds 


/ e l/3 c| + e * 7 si 


ml +2A ~ 


^ m LSND 


V 


0 e * 7 - e ifi 0 \ 

0 0 0 
• e * 7 c| + e l/3 0 

0 j 


(60) 


Opposite CP parities (a Pseudo-Dirac structure) of z / 3 and u 4 will again lead to (m u ) eT = 0. 
Now consider c v ~ s v . The mass matrix takes the form 


( e l/3 Cq + e* 7 . 




ml +2A 


r*j 



2 

LSND 


A 

0 


V 


(e * 7 - e ifi ) 

A 

(e 17 c| - 1 - e l/3 s|) 


(e * 7 - e i/3 ) \ 


e i 7 c| + e 


:7/3 Sq) 


( ( 

(e n c| + e 1 ' 3 s|) ) 


(61) 


We therefore find an approximate r-s symmetry, in analogy to the successful /i-r symmetry 
of the three neutrino case EDI- It is present when sterile neutrinos participate equally in 
solar and atmospheric neutrino oscillations. The r-s symmetry does strictly speaking only 
say that (m v ) TT = {m v ) s8 , here it holds in addition that (m u ) TS = {m v ) ss . If we consider 
the matrix 


m v 


( a b d d \ 

• h e e 

■ ■ / / 

■ • • f ) 


(62) 


we see that one eigenvalue is zero. Setting for simplicity b — h — e — 0 (these entries 
are suppressed in the previous equation), leads to two vanishing mass eigenvalues and 
777.34 = |(a + 2 / a/8 d? + (a — 2 /) 2 ), and therefore Ami = (a + 2 /) A/ 8 d 2 + (a — 2 /) 2 . 

In this limit the atmospheric Am 2 is vanishing. We have t/ e i = D e 2 = 0, |D e3 | J = | — 
(a/2 — /)/(-/ 8 d 2 + (a — 2/) 2 ) and I/^i = I/ r i| = 1/-/2. Hence, the approximate form of 
Eq. (f5Bj) is almost reproduced when a ~ —2/. Small breaking terms can in principle help 
to reach full agreement. 

The second interesting case (ii) occurs when m 4 ~ m 3 m 2 — m 4 ~ Am 4 ~ 0.1 eV. The 
implications are similar to case (i), but for completeness we give the resulting form of the 
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mass matrix: 


/ e l/3 Cq + e * 7 s 


ml +2A ~ 


^ m LSND 


O 


A c v (e * 7 — e l/3 ) 

^ £77 ^atm ^atm 

cl (e v 


c e + e 1 ' 


s e) 


V 


s,, (e * 7 — e l/3 ) \ 


Cr/Sr/ (e n Cq + e l/3 s|) 

s, 2 , (e i7 c| + eilS 4) ) 


.(63) 


" 5 ? ~© ' 0 °©7 

Comparing with case (i), we see that the second row of the mass matrix differs. It vanishes 
to first order when 77 = 0. If 77 7 ^ 0 the atmospheric neutrino mixing angle has some 
dependence on the form of the mass matrix. As an additional approximation, let us take 
9 q = 7 t/4 . Then, for e * 7 + e l/3 = 0, i.e., a Pseudo-Dirac structure of the two heavy masses, 
and c v — s v , we have 


*4 +2A ~ 


^ m LSND 


/ 0 

0 1 



0 1 

i 


■ 0 

0 

V ■ 


0 J 


(64) 


This matrix conserves L e -\- — L r — L s . Indeed, this global symmetry has been used in 

[B4i and, in somewhat different form in 165j . to explain the neutrino data including LSND. 
Moreover, the above matrix has all diagonal entries zero, a property typically shared by 
radiative models of neutrino mass generation. In Ref. |SBJ such a case is treated. 


C.2 The Mass Matrix in Scenario 2+2B 

For scenario 2+2B we can express the mixing matrix by exchanging in the mixing matrix 
from scheme 2 + 2 A the indices 1 3 and 2 4. Hence, 

/ cos 0q sin 9 e A A \ 

A A cos 6 > atm sin 6 > atm 

— cos 7] sin 9 q cos 77 cos 9 & — sin 77 sin 0 atm sin 77 cos 6 * atm 

y sin 77 sin 9 & — sin 77 cos 9 & — cos 77 sin 0 atm cos 77 cos 9 atm J 

Here we again put terms of order A ~ 0.1, which is the typical order of both the CHOOZ 
angle and the LSND parameter. We have again two cases of interest: 

(i) a/A m^ SND ~ 7774 - rn 3 > rn 2 ~ y/Ar 


U 


2+2B 


P . (65) 


vmf, » 777-1 , 

(ii) A / At77l SND — 7774 ~ ?77 3 7772 — 777! ~ 0.1 eV . 

Let us start with case (i), for which 

/ 0 A 


( 66 ) 


777 2+2B ~ 


^ m LSND 


Sr/ A 


Crj A 


J/3 


+ 4m e * 7 c atm s atm s v (e lf) - e n ) 


0 77 V'■'atm 


^ + 4m e i/3 ) 


CatmS at mC r? (e * /3 - e* 7 ) 
Cj}S,j(c. 


„2 

-"atm 1 


11 + 4m e-4 


c 2 f c 2 e 

V’-'atm ° 


3 atm c 

17 + 
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In contrast to scenario 2+2A it is the solar neutrino mixing angle whose contribution to 
the mass matrix is suppressed. Note again the approximate r-s exchange symmetry of the 
mass matrix in case of ,s v ~ c v . The (close-to-)maximality of 0 atm allows to further simplify 
the mass matrix to 


m 


2+2B 


\/ Am L 


LSND 


/o 

A 

S 77 A 

A 


e i/3 + g*7 

— e* 7 ) 

c^(e l/3 — e* 7 ) 



s 2 (e* 7 + e A ) 

c v s v (e 11 + e A ) 

V ■ 



c 2 (e* 7 + e l/3 ) 


(67) 


Opposite (identical) CP parities of u 3 and z / 4 lead to a vanishing /j/i entry and ts block 
(jit and /is elements). If we indeed impose a Pseudo-Dirac structure on z / 3 and z/ 4 , then 
we have 


/o 


m 


2+2B 


Am LSND 


V • 


A 

0 


S<q A C<q A \ 

5^ 


0 0 
■ 0 / 


( 68 ) 


If c v ~ s,, and the e/i entry of (recall that there can be a coefficient) is more suppressed 
than the er and es elements, then we have again a mass matrix conserving L e +L^ — L T —L s . 
Setting rj = 0 (atmospheric-sterile oscillations) and 0 atm = 7 t/4 , then at leading order 


/ o 


m?+ 2B ~ 


Am LSND 


0 

0 i(3 _|_ gi'T' 


0 0 \ 

0 e i/3 — e* 7 

0 0 


V • 


e i/3 + e i7 ) 


(69) 


conserving L e + L r . Recall that rj = 7 t /2 in scenario 2+2A lead to conservation of L M + L r . 
By choosing opposite or identical CP parities one can further simplify the mass matrix. If 
r) = 7t/2 (solar-sterile oscillations), then at leading order 


/o 


m 


2+2B 


Am LSND 


V • 


o o 

e i/3 _|_ e *7 gi/3 _ 

e ip + g®7 


°\ 

0 

0 

0 / 


(70) 


conserving L e + L s . 

Finally, we note that case (ii), defined by m 4 ~ m 3 S> m 2 — rn\ ~ Am 4 ~ 0.1 eV, leads to 
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